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The work presents the detailed analysis of the water dimer properties. Their parameters are 
investigated on the basis of a multipole interaction potential extended up to the quadrupole- 
quadrupole and dipole- octupole terms. All main equilibrium parameters of the dimer are ob- 
tained: its geometry, ground-state energy, dipole and quadrupole moments, vibration frequen- 
cies, etc. They are thoroughly compared with those obtained in quantum chemical calculations 
and from spectroscopic data. The efficiency of the present model potentials is discussed. A 
new viewpoint on the nature of the hydrogen bond is presented. The results of studies are 
thoroughly compared with the spectroscopic and computer simulation data. 

Keywords: water dimer, multipole interaction potential, dipole moment, quadrupole mo- 
ment, hydrogen bond 

moment Dd should be 6*o = 50° and Dd ~ 2.8 D, 
respectively. 

This brief review of the dimer properties shows 
that, nowadays, there is no consistent approach to the 
problem of formation of dimers, as well as to the cal- 
culation of their equilibrium parameters. The biggest 
difficulty is related to the absence of a clear approach 
to the magnitude of characteristic contributions to 
the interaction energy. 

In the present work, the water dimer properties are 
studied on the basis of a strictly defined potential that 
consists of: 1) literature-based dispersive and repul- 
sive interaction and 2) multipole electrostatic inter- 
action that is extrapolated to the overlapping region. 
Under the sufficient proximity of water molecules that 
leads to the overlap of the electron shells, it is nec- 
essary to use the quantum chemical calculations to 
determine the interaction energy. It was shown in [3] 
that the energy of a hydrogen bond itself does not ex- 
ceed ~ k^T m . The last circumstance is in agreement 
with the fact that the intramolecular distances in a 
water molecule change not more than by (1 ~ 3)% 
0, HI- This means that we can neglect the contri- 
butions of hydrogen bonds to the interaction energy 
of two water molecules. For the first time, this fact 
had been understood in the works of Berendsen [6] 
and Barnes 0. A critical review of the modern po- 
tentials, in which the existence of hydrogen bonds 
is ignored, is presented in Q. In the same work, the 



1. Introduction 

The main object of our research is the water dimer 
properties, which parameters are very sensitive to the 
type of intermolecular interaction. A water dimer 
is a closed system of two water molecules connected 
with the hydrogen bond. The equilibrium dimer con- 
figuration that corresponds to the minimum of the 
interaction energy is presented in Fig. 1. We will 
briefly discuss the main dimer parameters presented 
in Table 1 according to the quantum chemical calcu- 
lations and experimental data. In two types of 
dimers are studied: a "frozen" dimer, which has the 
distance between oxygen and hydrogen and the an- 
gle between O-H bonds fixed, and a "relaxed" dimer, 
where changes of the O-H length and the angle be- 
tween bonds are allowed. In [2], a "frozen" dimer con- 
figuration was used to determine the H-bond energy 
and the dipole moment of a dimer (the intermolecu- 
lar distance Too an d the angle 6o were postulated). 
From the comparison of the quantum chemical re- 
sults and the experimental ones, it follows that the 
optimum value of the ground-state energy of a water 
dimer is $ d » -(9 10)fc B T m , where T m = 273 K 
is the melting temperature of ice, and is observed at 
ro°c! = 2.98 A. The values of angle 6* and dimer dipole 
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Table 1. Characteristics of the equilibrium state of a water dimer according to different approaches 





Specification 


(0) s 


00, deg 


Ed 


Dd, D 
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J_ I (JZit.Il ttllllcl 


2.954 


19.3 


—9.94 




Relaxed dimer 


2.896 


20.5 


-12.78 






6-31G 


2.98 


60 


-10.69 


2.2 


m 


STO-3G 


2.98 


60 


-9.39 


1.72 


4-31G 


2.98 


60 


-14.18 


2.60 


w 




2.925 


51.8 


-9.11 




[10] 






-10.32 








nr Vj-ol 1 \j 


9 81 A 




— ±O.UD 


A 37 




HF G-311G++G(dp) 


3.001 




-8.85 


3.45 




HF G-311G++G(3df, 3pd) 


3.026 




-7.37 


2.88 


mi 


MP2 G-311G 


2.797 




-16.82 


4.24 


MP2 G-311G++G(dp) 


2.914 




-11.21 


3.30 




MP2 G-311G++G(3df, 3pd) 


2.904 




-9.79 


2.68 




BLYP G-311G++G(dp) 


2.927 




-10.03 


3.09 




BLYP G-311G++G(3df, 3pd) 


2.944 




-8.22 


2.54 


[12] 










2.6 


[13] 




2.94 






2.15 


® 


Exper. 


2.976±0.004 


57±10 








Exper. 






-9.96±0.4 




[14] 


Exper. 








2.60 



generalized Stillinger-David potential that eliminates 
its weaknesses and retains its positive features is pre- 
sented. Let us mention the attempt in Q to model 
the hydrogen bonds, but no results that deserve a sig- 
nificant attention are obtained. Specifically, the fol- 
lowing aspects are studied in the present work: 1) the 
ground-state energy of a water dimer along with its 
dipole moment as functions of the intermolecular dis- 
tance and the angles that define a relative orientation 
of molecules; 2) the influence of dimer's rotation on 
its ground-state energy and 3) vibration states of a 
water dimer. The possibility to use the electrostatic 
multipole potential for the description of attraction 
between water molecules in liquid water and the re- 
lations between the results obtained, computer simu- 
lations, and experimental data are discussed. 

2. Structure of the Interparticle Interaction 
Potential Between Water Molecules 

The interparticle interaction potential is modeled 
with the expression 

$(r, n) = $ r (r, n)+$ D (r, 0)+$ M (r, 0)+$ H (r, 0), (1) 

where $ r (r, f2) is the repulsive term, <&£>(r, O) de- 
scribes the dispersive forces, $m(»", ^) is a part of 
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the multipole expansion of the interaction energy be- 
tween two water molecules, and <&h(7, is the irre- 
ducible contribution caused by the overlap of the elec- 
tron shells of the water molecules (f2 is the set of an- 
gles that describe the orientation of water molecules). 
The reference to the multipole expansion for the in- 
terparticle interaction is supported by the following 
facts: 1) quantum chemical calculation of the multi- 
pole moments is a well-posed problem; 2) comparison 
of the different multipole contributions allows us to 
control the applicability of the multipole approxima- 
tion for the electrostatic interaction. Due to the wide 
discussion of the parameters of dimers 



17j . the last fact becomes stronger. 

In the present work, the multipole contributions 
to the interparticle potential are considered up to the 
three-moment terms, i.e. up to the octupole-octupole 
effects: 

<Mr)= E ^(Q^Q^'-D^-. (2) 
z — ' n\m\ r 

l<n,m<3 

Here, is the fc-th multipole moment, 

whose components are defined as Q^ k > =>• 

— — "1 "4 "ft 

is the differential operator and (:) stands for the 
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Fig. 1. Molecular coordinate system 




Fig. 

LCS 



2. Angles that define the orientation of molecules in 



sum over the indexes. The symbol (Q^ Q^' ') is 
used to determine the Kronecker product of the 
matrices that correspond to the multipole moments, 
and r is the distance between oxygens. After all the 
necessary differentiations, we substitute: r =>■ r OG 
and x a =>• r 00 ^2a, which corresponds to a selected 
dimer configuration. Multipole contributions include 
terms up to the dipole-octupole ones: 



$mW = <&DD + $ 



DQ 



QQ 



DO 



(3) 



Table 2. Components of the quadrupole 
moment of a water molecule in MCS 







Qyy, D A 


Qzz> D A 


[1] 


2.57 


0.37 


-2.94 


[3J 


2.34 


0.43 


-1.17 


[13] 


-0.09 


-1.67 


1.75 


[18] 


-0.07 


-1.61 


1.69 


[23] 


-0.19 


-1.58 


1.77 


[24] 


-0.07 


-1.61 


1.69 


[25] 


-0.55 


-0.84 


1.39 


[26] 


-0.07 


-1.74 


1.81 



which take the form 



r 3 
oo 



DQ 



z ' oo 



QQ 



oo 



$DO = -: 



z 'oo 



^ tt a W a/313 + U Q ^atf/V 



+ tfW<9 (2) + ^ 2 >0 (1) ) 

T a a u a22^ a a xy a22)- 



The components of the dipole moment of a water 
molecule in the laboratory coordinate system (LCS) 
are given by the relations 

d (1) = d(0, cos((5/2 + x), sin(£/2 + x )), 
d (2) = d(O,cos0,-sin0), 

where d is the absolute value of dipole moment of a 
water molecule, 5 is the angle in the molecular co- 
ordinate system (MCS) that defines the positions of 
hydrogens. It is stated that the water molecule is 
oriented in MCS as in Fig. 2. 

The values of components of the quadrupole and 
octupole moments, calculated in computer experi- 
ments or obtained experimentally are presented in 
Tables 2 and 3. The orders of magnitude for the mul- 
tipole moments are: 



d ~ 1 D, Q ~ 1(T 8 D • cm and O 



1CT 16 D-cm 2 . 



Here, D equals 1 Debye - unit of dipole moments. 
In our opinion, the optimal values of components of 
the quadrupole moment of a water molecule are ob- 
tained in the computer experiment [l3l | . These values 
are in good accordance with the experimental data 
(l8| . But the components of the quadrupole moment 
that correspond to the charge distribution in model 
potentials [1 , lj^- 22} differ even by the sign (see Ta- 
ble 2). The same situation is also characteristic of the 
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octupole moments (see Table 3). The mutual coher- 
ence of different works is only observed in the values 
of dipole moment. This fact gives us the ground to 
conclude that most of the model potentials are un- 
able to reproduce fine details of the intermolecular 
interaction. The transition between MCS and LCS is 
standard: 

d(O) = Rad, 

Qap($l) = RnQapRn, 
Oy a p(^) = RnOYapRji, 

where Rq indicates the rotation matrix by the an- 
gle fl: ft — a, 6, <p (a, 9, 4> are the angles of rotation 
around the LCS axes, as given in Fig. 3) 

For example, the rotation around the Oy axis is 
given by the direct and transposed matrices: 

(cos a sin a\ I cos a — sin cA 

o io, Rl=\ oi o • 
— sin a cos a/ \sina cos a / 

We note that the components of the quadrupole 
moment of a water molecule in [5, 15-20] were 
calculated with the use of the formula = 

5Ei(^«'4 ~ r i^ap) that differs from our one 
by the multiplier 3/2. 

According to the selection of LCS (see Fig. 1), 
we can use only Oy a p components of the octupole 
moment. By the same reason, as we use traceless 
quadrupole moments, we will use traceless octupole 
moments: 

OyIiP = Oy a p ~ -^OySap, 

where Oy — 53<v=i Oy a a is the trace of the matrix 
Oy a p- From our point of view, the most convenient 
values of Oy a p are obtained in [l3[ and where 
the values of components of the octupole moments 
are almost the same. 

The relative values of the multipole contributions 
of different orders to the intcrparticle potential are 
presented in Fig. 4. Here and below, we use the the 
dimensionless values for interaction energy $(r, fl) — > 
— > $(?•, f2)/fcBT m , where T m is the crystallization 
temperature for liquid water. The curves in Fig. 4 
correspond to the fixed value of angle Xt which is 



2 




Fig. 3. Role of the corresponding contributions to the in- 
terparticle potential: 1 - dipole-dipole interaction, 2 — dipole- 
quadrupole, 3 — quadrupole-quadrupole, 4 - dipole-octupole, 
5 — total multipole interaction 




-0.3' 4 3 

Fig. 4- The behavior of the sum of repulsive and disper- 
sive terms for the Buckingham (1), Buckingham-Corner (2), 
SPC (3), SPC/E (4), TIPS (5), and TIP3P (6) potentials 



supposed to equal 0.7° and does not depend on the 
distance between water molecules. The angle 9 is a 
function of the intermolecular distance and was ob- 
tained from the minimum of the interaction energy 
at each value of roo- The total value of multipole 
contribution is shown by curve 5. As we see, at a dis- 



Table 3. Components of the octupole 
moment of a water molecule 





O&x. D k2 


0^ YY , D A 2 


W XX , D A 2 


[11 


1.91 


0.30 


-2.22 


[13] 


-1.29 


-1.73 


3.03 


[18] 


-1.33 


-1.82 


3.16 
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Fig. 5. Profile of the potential surface <& = ^/k^Tm of the 
interaction of two water molecules ( T m is the temperature of 
melting for ice) 

tance of 3 A that is supposed to be the equilibrium 
for a water dimer, the depth of the total multipole 
interaction reaches —15. 

We should consider the behavior of the $do term 
that describes the dipole-octupole interaction. From 
Fig. 4, it follows that the curves of $do and &qq 
cross at r ~ r*, where r* » 3.2 A. This fact indi- 
cates the breaking of the multipole expansion. At the 
smaller distances, the electron shells begin to overlap, 
so it is the natural border for the usage of model elec- 
trostatic potentials. 

3. Repulsive and Dispersive Interaction 
between Water Molecules 

We will use the Buckingham-Corner (28j, three-point 
SPC [HIU, and TIPS [H H2 potentials to de- 
scribe the repulsion and dispersive effects between 
water molecules. In the Buckingham-Corner poten- 
tial, the effects of repulsion and dispersive forces are 
described with the sum of pair contributions of the 
hydrogen-hydrogen, hydrogen-oxygen, and oxygen- 
oxygen types: 



$(r) = Be- 3 ^ 



where B — (— e + (1 

As = /3r^i 6 , = 



'A 



Art' 



■-)e-*, A 6 = 



(4) 



p(l+/3)-8/3-6' 

e is the potential well 



depth, and r m is the intermolecular distance at the 
minima of energy. 

In the SPC and TIPS potentials, the dispersive and 
repulsive interactions inhere only in the electron shells 
of the oxygen atoms and are determined in the form 
similar to the Lennard- Jones potential: 



<J>(SPC) 



B, 



(SPC) 



-4, 



(5) 



The comparative behavior of the repulsive and dis- 
persive contributions to the Buckingham, Bucking- 
ham-Corner, SPC, TIPS, and TIP3P potentials is 
presented in Fig. 4. 

As we can see, the minimum of the interaction en- 
ergy is around 0.2 4- 0.3, and its position is approxi- 
mately at (3.5 -i- 3.6) A. Curve 1 that corresponds to 
the Buckingham potential differs very much from the 
latter potentials. Of course, the simple comparison of 
the dispersive and repulsive curves is not enough to 
select the appropriate potential. 

4. Ground State of a Dimer 

In this section, we will present the results of studies 
of the ground state of a water dimer, based on po- 
tential (1), where we will ignore the contribution of 
the short-range H-bond potential. We will consider 
the positions of the oxygens and hydrogens in water 
molecule remain still for the ease of calculations. Ac- 
cording to [Til], such requirement leads to the error 
not more than (1.5 3)%. The ground state of a 
water dimer is identified with the minimum of the in- 
teraction energy of two water molecules oriented, as it 
is shown in Fig. 1. More specifically, the equilibrium 
distance r °d between the oxygens and the angle 9q 
between the directions of unexcited dipole moments 
are found from the condition 



„(0) 



#o -H> absolute min $(f,6*,a = 0), 



where $ is the intermolecular interaction potential. 
It is considered to be a function of the dimensionless 
distance f — r 00 /v on , r OH — 0.97 A is the length of 
the O-H bond, and the angle a describes the rotation 
around the H-bond (see Fig. 2). 

The general view of the potential surface $(f, 6) is 
presented in Fig. 5. 

On the set of points of minimum for $(f, 8), the an- 
gle 9 m is a function of f. Its dependence is presented 
in Fig. 6. 
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The radial dependence of the interaction energy at 
the angle 9 m corresponding to the absolute minimum 
is presented in Fig. 7. 

It follows from Figs. 5-7 that the parameters pre- 
sented in Table 4 correspond to the absolute mini- 
mum of the interaction energy. 

One can see that the equilibrium parameters of a 
dimer are very sensitive to the selection of the form of 
repulsive and dispersive interactions. More precisely, 
these parameters depend on the law of decrease of 
the repulsive interaction. The depth of the potential 
well becomes bigger when the maximum of repulsion 
moves to the left (see Fig. 4). The position of the 
repulsive branch of the potential is a key factor that 
defines the equilibrium parameters of a water dimer. 
The value of the dipole moment of a water dimer is 
given by the formula 



D d = 2d w cos~(6 + {S/2- X ))- 



(6) 



The values given should be compared to the equilib- 
rium values obtained from the computer simulations 
and experiments (see Table 1). As follows from Table 
1, the intermolecular distance deviates from the lit- 
erature data. Another important dependence is the 
dependence on the angle a of rotation around the H- 
bond (see Fig. 9). The variable a is an intradimer 
characteristic, so it is responsive for the rotation of 
molecules inside the dimer. 

Another important characteristic of a dimer, which 
should be compared with other results Ill,ll2,ll7|, is 
the dependence on the angle 9 corresponding to the 
transversal vibrations of H-bonds (see Chapter 4). It 
is presented in Fig. 8. We see that, in the region 
of the minimum, the dependence of the interaction 
potential has a slope. As a consequence, in different 
models of water, the deviations of the equilibrium 
value of 9 can be huge. Moreover, the heat vibrations 
of the H-bond can influence the experimental values 
of the angle 9. 

It is necessary to mention that the most convenient 
dimer parameters to compare are: 1) the value of its 
dipole moment; 2) the values of components of the 
quadrupolc moment and it average value, which is 
defined as Q = \{Qxx + Qyy + Qzz) (here, we 
use the components of the non-traceless quadrupole 
moment), and 3) the vibration frequencies of water 
dimer. According to ([5]), the experimental value of 
the dimer dipole moment allows us to control the 




, . , . — - — ^ r ,A 

3 4 5 6 7 

Fig. 6. Radial dependence of the angle 9 m (5). Angles a, 0, 
and x remain fixed and correspond to the absolute minimum 
configuration. Curves 1—4 describe the dependence of 8 m in 
DD, DQ, QQ, and DO approximations, respectively 




Fig. 7. Dependence of the interaction energy on the inter- 
molecular distance in a dimer configuration 

Table 4- Equilibrium parameters of a water 
dimer for different interaction potentials 







Ed, ksTm 


6 


X 


D d , D 


B 


2.96 


-8.5 


28.03 


0.75 


2.8435 


BC 


2.89 


-9.65 


28.70 





2.8148 


SPC 


2.72 


-12.49 


27.69 


1.03 


2.8143 


SPC/E 


2.72 


-12.49 


27.69 


1.03 


2.8143 


TIPS 


2.7 


-12.66 


27.61 


1.14 


2.8137 


TIP3P 


2.69 


-12.90 


27.60 


1.16 


2.8136 



equilibrium values of the angle 9. The components 
of the quadrupole moment of a water dimer are con- 
nected with the components of the quadrupole mo- 
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Fig. 8. Comparative behavior of the interaction potentials 
(1) - B, (2) - BC, (3) - SPC, {4) - TIPS 



-l<r 

Fig. 11. Dependence of the interaction energy on the angle \ 




Fig. 9. Dependence of the dimer energy on the angle of rela- 
tive rotation around the H-bond 




-10 



Fig. 10. Dependence of the interaction energy on the angle 
at t-oo = 2.9757 A 




Fig. 12. Temperature dependence of the equilibrium distance 
between oxygen atoms in a water dimer 



O' 

Fig. 13. Vibration modes of an H-bond 

merit of an isolated water molecule by the relation 

= C?S + «S + l(R*(4 w} - d[ w \ + 

+ Rp^-dM) a ), (7) 
The average value of quadrupole 



where R = r, 
moment is 



(0) 
oo ■ 



1 



)d = Qi+Qi + ^dto^oo (cos 6 - cos 5/2). 



(8) 



It is worth to mention that the components of the 
dipole and quadrupole moments can be easily cal- 
culated with the help of the methods of quantum 
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Table 5. Frequencies of small oscillations 





B 


BC 


SPC 


TIPS 


"Frozen dimer" 


Rotating dimer 


"Frozen dimer" 


Rotating dimer 




cm" 


l 


175.12 




166.23 


182.51 




164.20 


63.14 


68.91 


1^2, 


cm" 


l 


196.82 




194.44 


198.14 




194.43 


244.63 


212.63 


<*>3, 


cm" 


l 


391.83 




385.63 


400.87 




386.73 


246.00 


247.64 


Co>4, 


cm" 


l 


521.64 




513.29 


529.97 




511.42 


307.93 


313.00 



chemistry. Another surprising fact is the possibility 
of the existence of metastable states of water dimers. 
With the fixed angle 9 that corresponds to a min- 
imum, there is a possibility of the turn of another 
water molecule on a specific angle, that can be cal- 
culated studying the dependence of the interaction 
energy on angle \- 

As follows from Fig. 11, there are two local min- 
ima separated by the potential barrier. The depth 
of the second minimum is approximately 5k B T m , and 
its position corresponds to the reflection of a molecule 
relative to the xy plane. 

5. Rotation of a Dimer 

If a dimer is formed in the gas phase, its thermal 
motion consists of the translational, rotational, and 
vibration modes. It is clear that only rotations and 
vibrations influence the equilibrium parameters of a 
dimer. Let us mention that the experimental val- 
ues of intermolecular distance @,[l,[l3| correspond to 
the rotating dimers. This circumstance is ignored in 
the computer experiments. In this section, we will 
consider only the rotational influence on the param- 
eters of a dimer. The rotation of a dimer is followed 
by the emergence of centrifugal forces, which lead to 
the growth of the equilibrium distance between oxy- 
gens. Rotation takes place around two axes that are 
perpendicular to the OX axis. Both axes lie in the 
molecular planes. The growth of the distance between 
oxygens is determined by the equation 

k{r - r ) = m (uf + uf)r, (9) 

where the average frequencies of rotation can be esti- 
mated from the relations 

2 k B T 2 k B T 

w l = -T-, w 2 = —j— ■ 



The inertia moments I\ and li are equal to: I\ — 
= 1.332 x 1(T 33 g • cm 2 , I 2 = 1.219 x 1(T 33 g • cm 2 . 
The values of force constant k are calculated in the 
next chapter. But, regardless the calculations, the 
equilibrium distance grows linearly with the temper- 
ature. The rise of the equilibrium distance between 
oxygens is presented in Fig. 12. 

As we can see, the distance between oxygens at 
the melting temperature increases approximately by 
0.04 A: from 2.96 A to 3.00 A. The ground-state en- 
ergy decreases by 0.04. In a vicinity of the critical 
point, the equilibrium distance reaches 3.06 A. 

6. Small Oscillations of a Dimer 

From the configuration of a dimer, it follows that 
the most characteristic types of small oscillations are: 

1) longitudinal one connected with a change of the 

r OQ distance between the oxygens: r ~ r °° ; 

2) two transversal oscillations of the H-bond that are 
connected with small rotations around the x (by angle 
9) and z (by angle (j) ) axes and 3) intradimer oscil- 
lations connected with a change of the relative orien- 
tation of water molecules during the rotation around 
H-bonds. Here, it is necessary to mention that the 
variables f and 9 are not independent and form a 
new pair of normal (hybrid) coordinates. Longitudi- 
nal and two transversal vibrations can also be inter- 
preted as the oscillations of the H-bond (Fig. 13). 
We have to remember that dimers rotate; therefore, 
the spectra of small oscillations will change. 

According to this, the Lagrange function for small 
oscillations of the dimer takes the form 

L = - TUijiiXj — - KijXiXj, (10) 

1<«J<4 1<«J<4 

where Xi = (f, 9, </>, a) are the generalized coordi- 
nates, and ii — (f, 9, <f> 1 a) are the generalized ve- 
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locities. It is easy to see that the mass tensor has a 
diagonal structure: 

mn = I r , m 22 = I e , m 33 = 1$, 77144 = 

where I r = /ir^. /i = l/2mw is the reduced mass of 

a water molecule, / = j^j—^rj are the components of 

the reduced inertia moments, i — 8, a, and the super- 
scripts number water molecules. In accordance with 
Fig. 1, the components of the inertia moment are: 

— o ra r 2 ■ 
lf ] = 2m H r 2 H cos 2 5/2; 
jW = m H r 2 H (l-|-cos 2 <5); 
I {2) = 2m H r2 H (l - cos 2 8/2 cos 2 6) 

and are equal to 
I T w 14.14 x 10~ 46 g • cm, 
Ig w 0.83 x 10" 46 g • cm, 
1$ w 0.61 x 10~ 46 g • cm, 
I a » 0.95 x 10~ 46 g • cm. 

The force constants k^ are determined in a standard 
way: h 



d A <i> 



v.jj — Qx-dx Ui=o- ^ u ^he derivatives are cal- 
culated at a fixed value of angle x- The frequencies 
of small oscillations are calculated in a standard way 
and are presented in Table 5 for the potential with 
repulsive and dispersive parts from the Buckingham 
and Buckingham-Corner potentials. Two first fre- 
quencies correspond to hybrid oscillations of the f, 6- 
type. Their normal coordinates are (at x = 0): 



Ui 



(fcn — Jj-wf )f + k 



12 1 



I r (uj 2 - u\)(k 22 - h^l) 



C\ exp(zwit), 



k 12 r + (k 22 - I 6 uj 2 )0 
u 2 = — -2 7 — 2\G 2 exp(m> 2 i). 



(11) 



(12) 



Coefficients Ci and C 2 are determined from the initial 
conditions. Thus, we can only speak about a single 
transversal oscillation of an H-bond, while the other 
one mixes with its longitudinal oscillation. 



7. Discussion 

From the qualitative point of view, the formation 
of dimers is related to the formation of a hydro- 
gen bond between two water molecules. In this 
case, the ground-state energy of a dimer is identi- 
fied with the energy of the H-bond. In our dimen- 
sionlcss units, it should be equal to -10. This esti- 
mate is confirmed by the quantum chemical calcu- 
lations 0, [l(| [l6[ . The close value of bonding en- 
ergy (Ed = —9.19) is obtained, by using a gener- 
alized Stillinger-David potential [16fl. However, the 
different approach is considered in the present work. 
It is supposed that the dimers are formed due to 
three well-defined types of interaction: repulsive, dis- 
persive, and multipole electrostatic interactions. In 
our calculation, we restrict ourselves to all contribu- 
tions up to the dipole-octupole ones. The extrapola- 
tion of the interaction to the overlapping region does 
not lead to any anomalous increase of the dipole- 
octupole interaction relative to other terms. This 
is the evidence of the legitimacy of using the de- 
scribed procedure of extrapolation on the distances 
up to r 00 = 3 A. The reliability of the results is 
confirmed by the fact that the values of quadrupole 
and octupole moments calculated in the computer ex- 
periments coincide with the experimental data. The 
optimal value of ground-state energy of a stationary 
dimer according to our calculations is E c i — —9.65 
and is reached at the distance ro°o = 2.89 A. It is 
shown that the rotation of a dimer, which is quite 
natural in the gas phase, leads to an increase of the 
equilibrium distance by 0.04 A. One more circum- 
stance that can influence the experimental value of 
distance between oxygens is the excitation of longi- 
tudinal oscillations that correspond to the frequency 
W\. By the order of magnitude, the amplitude of os- 
cillations of dimers is 



I At 



lk B T m 



0.1. 



We see that the rotation and oscillations can change 
the distance between oxygens by 0.1 A. This fact can- 
not be ignored. Another parameter of a dimer that 
can be controlled is the angle 6q, which defines the 
equilibrium mutual orientation of the dipole moments 
of water molecules that form the dimer. This angle 
is directly connected with the values of dipole and 
quadrupole moments of a dimer, which can be cal- 



286 



ISSN 2071-0186. Ukr. J. Phys. 2013. Vol. 58, No. 3 



Dimerization of Water Molecules 



culated with the help of quantum chemical methods 
or be obtained experimentally. In our work, the fol- 
lowing values of the angle and dipole moment of a 
dimer correspond to the above interaction energy and 
the equilibrium distance: 6 = 28.7°, D d = 2.81 D. 
These values slightly differ from those obtained exper- 
imentally or calculated within the quantum chemical 
methods (see Table 1). In [ll, 17 1, it was shown that 
quantum chemical calculations result in the dipole 
moment of a dimer Dd = 2.6 D. According to ([6]), 
this value corresponds to the angle 9q m 36°. The 
value of 6*o was also studied in [7, 8]. It was shown 
that, at r OQ = 2.976 A, the angle O equals 57° ± 10°. 
From the spectroscopic experiments, it follows that 
the angle should be 51° ± 20°. Such a difference in 
the values of 8 is related to the sloping character of 
the dependence in a vicinity of the minimum, as 
it follows from Fig. 8. Vibrations and the rotation of 
a dimer also influence the value of 9q: 



\A9\ 



2.6°. 



The influence of the interaction of water molecules 
on the parameters of water molecules will affect, in 
some way, the equilibrium parameters of a dimer. 
It was demonstrated in [5] that the dipole moment 
of a water molecule increases by 1 1.5%, when 
molecules approach it to a distance of 3 A. So the 
growth of the dipole-dipole interaction will result in 
changes of the intermolecular distance and the in- 
teraction energy. The overlapping effects that man- 
ifest themselves at the distances smaller than 3 A 
should not be considered because the relative shift 
of valence vibrations of a proton at the condensa- 
tion does not exceed 1 3% [14]. But the exis- 
tence of hydrogen bonds should not be neglected be- 
cause they influence the value of heat capacity (29| . 
The weak overlapping of the electron shells shows it- 
self in the tunneling of protons from one molecule 
to the other one along the hydrogen bond. But 
these effects are characteristic only at the superlow 
frequencies (~30 cm -1 ). All these qualitative ar- 
guments allow us to state that the energy of the 
hydrogen bond itself does not exceed k-QT m . The 
relatively small contribution of hydrogen bonds to 
the interaction potential was stated in the works by 
Dolgushin 3(3, 3l| . He showed that the sharing in- 
dices of the electron density between water molecules 
under the influence of a neighbor molecule are less 



than 3%. These results were confirmed by the later 
work by Fulton [H]. We mention that, in 0,0], 
hydrogen bonds were ignored without any justifica- 
tion. In those works, it was shown that the in- 
teraction energy calculated with the potential con- 
sisting only of the classical electrostatic and repul- 
sive interactions matches the quantum chemical re- 
sult with sufficient accuracy. The multipole approx- 
imation allows us to easily build the averaged inter- 
action potential between water molecules. It should 
be noted that certain limitations of the rigid multi- 
pole interaction will manifest themselves in the de- 
scription of multimer properties, where multiparti- 
cal effects take place. In particular, the equilibrium 
distance between oxygens reduces to r OQ — 2.78 A 
in liquid water. The denial of the model of hy- 
drogen bond as the specific donor-acceptor type of 
interaction finds its confirmation in another fact. 
The values of self-diffusion coefficients and shear vis- 
cosity for water have the same order of magnitude 
as those for the liquids that do not form hydro- 
gen bonds. This fact indicates that the character 
of translational and rotational motions in these liq- 
uids is similar. It cannot be agreed with the exis- 
tence of sharply directed hydrogen bonds with the 
interaction energy of ^\Qk^T m . Our approach does 
not have such complexity. The electrostatic multi- 
pole interaction satisfies the superposition principle, 
so the total electric field formed by the randomly po- 
sitioned and orientated water molecules is relatively 
small and does not affect critically the movement of 
water molecules. 
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II.B. MaxAaunyK, M.Y1. MaAOMywc, LB. TKmaHKiK 
flHMEPISAUM MOJIEKYJ1 BQUH. 

MOflEJIPOBAHHil MDKMOJIEKyjIflPHOi BSASMOflll 
HA OCHOBI MYJIBTHnOJlBHOrO ilOTEHUIAJiy 

P e 3 K3 m e 

Po6oTa npHCBHHyeTtCH fleTajibHOMy aHajii3y BjracTHBOCTeft 
flHMepiB BOflH. Bci niiTaHHH, noB'H3am 3 uieio npo6jie- 
mok>, flOcjiifl>KyiOTbCH Ha ocHOBi MyjibTHnojibHoro no-TeHuia- 
jiy, no,z],OB:>KeHoro KBaflpynojib-KBaflpynojibHoro Ta flunojib- 

OKTynOJIbHOrO BHeCKiB. OTpHMaHO 3HclMCHHH BCix OCHOBHHX 

piBHOBa:>KHiix napaivieTpiB .znnviepy bo^h: Horo reoMeTpiiMHHx 
xapaKTepiiCTHK, eHepri'i ochobhcto CTaHy, .zpincnbHoro i KBa- 
flpynojibHoro MOMeHi-iB, nacTOT KOjiHBajibHHx CTamB tohj,o. 
/leTajibHO o6roBopioeTbCH afleKBamicTb MOflejibHHx noTeimia- 
jiiB MijKMOjieKVjmpHOi B3aeMO^,ii" y BO^i. IlponoHyeTbCsi hobiih 
norjiflfl Ha npnpofly BO^HeBOro 3b'h3kv. Pe3yjibTaTH ^ocjiifl^Ke- 

HHH peTejIbHO 31CTaBJieHO 3 flaHIIMH CneKTpOCKOnlHHHX flOCJli- 

AJKeHb i KOMn ' K)TepHHx CHMyjiHH,iii. 
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